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Abstract —The physical-layer security issue in the multiple 
non-regenerative wireless-powered relay (WPR) networks is in¬ 
vestigated in this letter, where the idle relay is treated as a 
potential eavesdropper. To guarantee secure communication, the 
destination-based artificial noise is sent to degrade the receptions 
of eavesdroppers, and it also becomes a new source of energy 
powering relays to forward the information with power splitting 
(PS) technique. We propose an efficient algorithm ground on 
block-wise penalty function method to jointly optimize PS ratio 
and beamforming to maximize the secrecy rate. Despite the 
nonconvexity of the considered problem, the proposed algorithm 
is numerically efficient and is proved to converge to the local op¬ 
timal solution. Simnlation results demonstrate that the proposed 
algorithm outperforms the benchmark method. 

Index Terms —Energy harvesting, non-regenerative relaying, 
power splitting (PS), beamforming, multiple wireless-powered 
relays, physical-layer security, artificial noise. 

I. Introduction 

Simultaneous wireless information and power transfer 
(SWIFT) is a promising solution for prolonging the lifetime 
of energy-constrained wireless nodes [1]. Due to its efficient 
utilization of wireless spectrum to harvest energy and process 
information from the same wireless signal, SWIPT triggers 
a great deal of research interests in cooperative transmission 
[2]-[5]. Two protocols, namely time switching (TS) and power 
splitting (PS), were proposed for energy harvesting relay with 
SWIPT in [2]. TS enabled relay utilize different time blocks 
to realize information processing (IP) and energy harvesting 
(EH), while relay proportionally split the received signal power 
to accomplish these two operations separately in PS protocol. 
In [3], the authors investigated a multiple-antenna relay system 
with SWIPT where a “harvest-and-forward” strategy was 
proposed to maximize the achievable rate. Meanwhile, the 
authors in [4] investigated power allocation strategies and their 
impact on the multiple source-destination pair networks with 
an EH relay. And the performance of a large-scale network 
with multiple transmitter-receiver pairs with/without relaying 
was discussed in [5], where receivers employed PS technique 
to harvest energy. 
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Recendy, with the rapidly developed SWIPT technologies, 
physical (PHY)-layer security has attracted significant atten¬ 
tion in SWIPT with application to wireless networks. As is 
mentioned in [6], [7], EH receivers are generally deployed 
sufficiently close to the transmitter than information receivers 
(IRs) to achieve higher signal power. This gives rise to 
a challenging PHY-layer security issue that EH receivers 
may easily eavesdrop the confidential information with better 
channels, rather than harvest energy as they are presumed to 
do. To resolve this issue, substantial research efforts have 
been dedicated to SWIPT with the consideration of PHY- 
layer security [6]-[9]. The authors in [6] investigated secrecy 
wireless information and power transfer with joint informa¬ 
tion and energy beamforming design in multi-input single¬ 
output (MISO) system, while the dual use of artificial noise 
(AN) and energy signals for facilitating simultaneous secure 
communication and energy transfer with MISO beamforming 
was introduced in [7]. The authors in [8], [9] split the power 
at the transmitter to send confidential message to IR and 
an AN to interfere ER against eavesdropping. Moreover, the 
secrecy problem in cooperative networks with SWIPT was 
considered in [10]-[12]. The authors in [10] focused on secure 
relay beamforming for SWIPT in one-way relay systems with 
an external eavesdropper and EH receiver. Meanwhile, the 
authors in [11], [12] investigated a wireless-powered jammer 
system, where the jammer (e.g. idle EH receiver [12]) could 
harvest energy from the wireless signal and use it to interfere 
with the external eavesdropper. However, these works [10]- 
[12] associated with cooperative transmission did not consider 
that EH receivers may intercept the confidential information 
intended for IRs. 

Motivated by the aforementioned problem, in this letter, we 
consider a multiple non-regenerative wireless-powered relay 
(WPR) networks, which is different from the scenario in [3] 
with only one WPR. In this networks, WPRs are supposed 
to harvest energy to assist the source-destination transmission 
when active. However, when the relay is idle, it may become 
a potential eavesdropper in the network, and maximize its 
decoding signal-to-noise ratio (SNR), rather than split power 
for cooperative transmission. To improve the security of the 
source-destination transmission, the destination-based AN is 
applied to our scenario to degrade the receptions of idle 
relays, which also becomes a new source of energy powering 
WPRs in the second-hop transmission. Distinguished from the 
exhaustive search for PS ratio in [3], we propose an algorithm 
based on block-wise penalty function method to find the local 
optimal PS ratio and optimize beamforming for each WPR 
to maximize the secrecy rate with much less computational 
time. Despite the nonconvexity of the considered problem, the 
proposed algorithm is numerically efficient and is proved to 
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converge to the local optimal solution. Numerical experiments 
manifest the superior performance of the proposed scheme as 
compared with the benchmark method. 

Notation: We adopt the notation of using boldface for 
column vectors (lower case), and matrices (upper case). The 
hermitian transpose is denoted by the symbol (•)^. For a 
complex scalar x, its complex conjugate is denoted by x*. 
E[-] and CAf{-) denote the statistical expectation and complex 
Gaussian distributions, respectively. For vector x, diag(a:) 
means putting x on the main diagonal. 


scheme, pi denotes the ratio of power split for EH, and the 
left power for IP is {1 — pi) at r^. Hence, the parts of received 
signal for EH and IP at are written as follows, respectively 

yf^ = \/lh{hsiXs + hdiXd + ni), ( 2 ) 

Pi — \/1 Pi{hsiXs T hdiXd T 'Oi) rTcz, (3) 

where Ud £ C is the additive noise introduced by signal 
conversion from the passband to baseband with the distribution 
CA/'(0,(T^j) at Xi. The harvested energy of I'i is given by 


II. System model 

Consider a two-hop parallel relay networks where a source 
S, a destination V and a set of relay nodes {ri,i = 1, • • • , N) 
are all equipped with a single antenna. Each Xi can assist 
the source-destination transmission when active. In this paper, 
we assume that the active relay Xi can not only process the 
received signal, but also harvest energy from it to help S 
forward the information resorting to PS technique. However, 
when the relay is idle, it may become a potential eavesdropper 
in the network, and maximize its decoding SNR, rather than 
split the power for cooperative transmission. What’s worse, if 
the idle relay is closer to S or other active relays, it can easily 
eavesdrop the conhdential information with better channels. 

To guarantee secure communication, at the time of con¬ 
hdential signal transmission, 1) concurrently sends AN to 
the area of relay nodes to degrade the channel conditions of 
idle relays. At the second hop, Xi amplihes and forwards the 
received signal to V depending on its harvested energy. In 
this circumstances, AN can not only be used to degrade the 
receptions of eavesdroppers, but also become a new source 
of energy powering the information delivering from relays 
to 1) [7]. This scenario is potentially applicative for sensor 
networks where a sink node with conventional supply of power 
wants to send its data to sensors far away from it. So the sink 
node asks for the help of some intermediary sensor nodes. 
However, the intermediary sensor is not willing to help the 
sink node due to the limited battery capacity. Hence, SWIPT 
is employed to encourage those intermediary sensors to use the 
harvested energy to participate in cooperative communication 
[13]. Nevertheless, in this application, some idle sensors may 
prefer to eavesdrop the legitimate information, rather than help 
the cooperative transmission, which incurs serious challenges 
to the PHY-layer security. 

Denote the channel responses from S to Xi, from V to 
the active Xi and from Xi to V by hgi £ C, hdi £ C and 
hid £ C, Vi. At the first hop, S sends the legitimate signal Xg 
to relay nodes, meanwhile, V transmits AN Xd to confound 
the potential eavesdroppers in the area of relay nodes. The 
received signal at relay Xi can be expressed as 

Vi — hgiXg hdiXd T Xli^ (1) 

where £ C represents the additive white Gaussian noise 
(AWGN) associated with relay node Xi following the distri¬ 
bution CM{0,af), meanwhile, Xg and Xd are the transmitted 
symbols of S and 27 with the conventional energy supplies, 
such as E[xsX*] = Pg and E\xdX*^ = Pd- Based on PS 


Pf^ — VPi {Ps\hgi\'^ -f Pd\hdi\^ + (Xi) , (4) 


where p £ (0,1) denotes the energy conversion efficiency from 
signal power to circuit power. 

In order to get more insight of this relay networks, we 
assume that there is only one idle relay (r^v+i) eavesdropping 
the legitimate information in this letter, and denote it as 
the potential eavesdropper {£). The channels from S to £, 
from V to £ and from Xi to £ are given by hge, hde 
and hie,yi- Therefore, the received signal-to-interference-and- 
noise ratio (SINR) of £ at the first hop can be written as 
SINRe 1 = 2 , and is the AWGN at £. 

’ ^d\^de\ ~r^e ^ 

We denote Wi as the scalar beamformer at relay Xi, and 
utilize it to deal with the received signal at the corresponding 
relay node. Assuming perfect self-interference cancellation, 
the received signal of 27 at the second hop can be given by 

E n - 

y 1 PihidWihgiXg -f 

1—1 

E n - 

. ^\J^-PihidWin^py hidWpndPnd,{'£) 
2 = 1 ^ -^2=1 


where Ud is the AGWN introduced by the receiver antenna at 
27 with the distribution cr^). Eor notational convenience, 

we make of = a^i = a1 = , Mi in the following parts. 

Making use of matrix-vector expressions, the achievable rate 
of 27 can be written as 


Rd{w,p)=\ogJl + 


B pGiB pW 


BpQ^BpW + w^QiW + a'‘ 


, ( 6 ) 


where Gi = Pggig\, and g^ = [hidhgr,--- ]hNdhgN] £ 
C^xi. w = [wi;-- - ; wAf] £ p = [ pi]-■ ■ ; pat] £ 

Bp = diag(Vl-pi,--- ,y/l- pn), and = 
diag • • ■ , cr^- 

To maximize the eavesdropping rate, £ combines the two 
received signals in both stages using maximal ratio combining 
(MRC), and its achievable rate is given by 


2?e(tU,p) = log2 1 -f 


Ps\hg 


Pd\h 


de 


B pG2B pW 


w'i B pG^B piu+w'i B PQ 2 B pW+w'i Q2W-\-a- 


(7) 


where G 2 = Psg292 ^^d ^2 = [hiehgi,--- \hNehgN] £ 
Civxi; ^ PdPi^gl 93 = [fiiefidi; • • • ; 
and Q 2 = diag (|hiep, • • • , \hNe?) 
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III. Problem Formulation For Secrecy Rate 
Maximization 

In this section, we investigate the joint PS and beamforming 
design to maximize the secrecy rate Rgr, which is given by 

Rsr{w, P) = ^ [Rd{w, p) - Re{w, p)]+ , (8) 

where [a]+ = max(0,a), and i is from the fact that the 
source-destination transmission takes place in two time slots. 

For the considered scenario in this paper, the optimization 
problem can be written as 


max 

w.p 

s.t. 


Rsr{w,p) (9) 

[(1 ~ Pi)Ps\hsi\'^ -f (1 — Pi)Pd\hdi\'^-\- 

{ 2 -p^y] e [o,i],Vi = i,--- ,iv, 


where the transmit power of Vi is constrained by the harvested 
energy at ri, Vi. It is noted that the considered optimization 
problem (9) is non-convex, and the combined inequality con¬ 
straints further increase the difficulty to globally solve (9) 
with respect to the two coupled blocks of variables [w, p). In 
order to handle the nonconvexity and inequality constraints, we 
introduce the penalty function method [14] into our problem. 
Before that manipulation, formulation (9) should be equiva¬ 
lently transformed into the following problem with the newly 
introduced slack variable ip. 


max Rsr{w,p) (10) 

w,p,%p 

s.t. \wi\^ [(1 — Pi)Ps\h8i'\^ -f (1 — Pi)Pd\hdi'\^-\- 

{2-p,)a^] = e [0,1], 

lA, >0,Vi = l,--- ,N. 


Further let us dehne the following penalty function 

f{w,p,xp,X)=Rsriw,p)-X^ [{1-Pi)Ps\h8^f 

-f(l - Pi)Pd\hdi\'^ + (2 - Pi)cr‘^] - P™\^ , (11) 

where the second term represents the penalty term constituted 
by a total of N constraints in (10), A denotes the penalty 
parameter. Notice that problem (10) is equivalently rewritten 
as the maximization problem of (11) constrained by the 
feasible set of (10), in which the penalty term equals zero. 
If the obtained solutions are not feasible, which violate the 
constraints of (10), the penalty term may increase much with 
the quadratic growth, resulting in lower value of (11). The 
penalty function is introduced to balance the secrecy-rate 
performance against the penalty term through the setup of 
A. Therefore, the key idea of penalty function method is to 
increase the penalty parameter A progressively in order to 
obtain an optimal feasible solution of (10). 

In consideration of the problem structure, we further in¬ 
troduce Block Coordinate Gradient Descent (BCGD) Method 
[15] into our problem, which can be utilized to handle the 
large-scale multiple-block problem by computing and unifying 
each block variable through sequential updating. And the 
target of each update in our proposed algorithm is to increase 
the value of (11) by reducing the penalty term or enhancing 


TABLE I 

The Block-wise Penalty Function Method 


Initialization: Generate feasible block variables (w^, , 1 )}^) 

and penalty parameter A®. 

SI: For k = 1, - ■ ■ ,K, do S2-S3 until converge 
S2: Fori = I,-- - do 

(1) = w'p + ,\*‘), Vi; 

(2) = pfc -P/3p.<^p.(^ofc+^pfc,),-0^A'=), Vpi e [0,1]; 

(3) >0. 

Note that and 7 ^. are the step sizes in each update, 

and can be obtained by 

Armijo stepsize rule: For the considered maximization problem, 
choose 0 < < 1 and G (0,1). Let be the largest 

element in ^ satisfying 

/(3^fe+i) _ /(ajfe) > 

where can be replaced by w^, and , Vi, and V/(ai^) 
turns to be the partial derivative of / with respect to i.e., 
^(x^). And is the ascent direction for our problem. 

. is set as follows 


The ascent direction d^u 


(m^(, pV A'") = -(mf( , r/,^ A'=). 

Similarly, the other two , p'j^yy , \^) and 

d^.{'w^+^, jipyy, can be achieved. 

S3: Stopping Criteria (Sequential Convergence) 

IIP'=+MF 


If min 


„fc||2 ||pfc + l_pfc||2 ||^fc + l_.^fc||2 

‘ ||r/)'' + l||2 


jpuJT+Tp" 

Eil'i [{l-pt")Ps\h 8 i\^+{l-P^+^)Pd\hd^\^ 

< e, 


-P(2 - p'^+^W] + 


^ k+l _ pEH 


Stop and return (io''+V p'^+V •i/;''+l); 

Else Set A*^+^ = cA*(c > 1), fc := fc + 1 and go to SI. 

End 


Rsr- With the help of Armijo stepsize rule (refer to Table I), 
the monotonic increment of (11) is guaranteed in each update. 
When the stopping criteria is met according to Table I, the 
constraints of (10) are all satished, making the penalty term 
equal zero, and the local optimal solutions are obtained. 

The block-wise penalty function method is a combination 
of penalty function method and BCGD method, the details of 
which are given in Table I, where e is the tolerant error, 
p^yy and '4’\i) 3IS the temporary vectors with latest updated 
information instead of p^ and 



= K+^-' 


■ ■; 'd’%) 

Pm 

= 

'iPi—liPi') Pi -\-1 ’ * * * 

; Fa)) 


= (01+'; •• 

•;0fj^;0f;0Li;-- 

■ ;07v)- 


Before presenting the convergence result, we first assume 
the feasible set is nonempty, which means that there is at least 
one feasible point of (10). Next we have the following result; 

Theorem 3.1: Any limit point of the sequence 
'0^} generated by the proposed block-wise penalty 
function method is a Karush-Kuhn-Tucker (KKT) point of 
problem (10). 

Proof. First let (tu*,p*,'0*) be any limit point of the se¬ 
quence p^', 0*^}. Due to the property of Armijo stepsize 
rule, we have the monotonic increment of /(lu^, p^, 0^, A^), 






4 


and obtain the following inequality 

= Rsr{w\P>^ - p\ -0'=). (12) 

Rearranging (12), we have 

, p'^-fiw^, p°, 0°, A°) <00,(13) 

where p(w^, p^,il)^) represents the second term without A in 
(11). Considering the boundedness of Rsr and the property of 
the limit point {w*, p*, 0*), we have X°°p{w*, p*, 0*) < oo. 
Because X°° —?► oo, we can deduce p{w*, p* j'ljj*) = 0 and 
(lu*, p*, 0*) is a feasible point of (10). If {w*, p*, 0*) is not 
a KKT point, there should be another feasible point {w, p, 0) 
making 

Rsr{w,p)> Rsr{w*,p*). (14) 

Further, let us consider the subsequence { ki }, where 
lim {w ^', , 0^') = (lu*, p*, 0*), so that 

i—¥oo 

lim f(w'^' ,p^' ,0'=* ,A) =f(w*,p*,ip*,X)=R,r(w*,p*),yX, 

where (a) is from the fact that the penalty term of (11) equals 
zero at point {w*, p *With a view to (14), for large 
enough ki , we have 

f{w, p, 0, X'^i=Rsr{w, p)>Rsr(:w*, p*)>f(}V^%p^%xl)^\X^'), 

which is a contradiction to the fc^-th iteration with respect to 
A^% where (b) is due to the fact that {w,p,ip) is a feasible 
point of (10) making the second term of (11) equal zero. □ 

IV. Simulation Results 

In this section, we investigate the performance of the 
proposed algorithm (PA) via numerical simulations. The path 
loss model for the energy harvesting relay channel is Rayleigh 
distributed and denoted by where |/3| and d represent 

the short-term channel fading and the distance between two 
nodes. |/3p follows the exponential distribution with unit mean 
[3]. The position of each node is given in Fig.l. We set 
the noise power as = —40dBm, the energy conversion 
efficiency p = 0.5, the penalty parameter A = 35, the constant 
c = 1.2 and the initial PS ratio pi = 0.1, Vi for all simulations 
results, which are averaged over 100 channel realizations. 

A benchmark method is set as a comparison, named as 
simple amplify-and-forward (SAF) scheme, where we per¬ 
form exhaustive search to obtain the power splitting ratio 
for each relay, and apply a simple amplified scalar 6 i with 
the consideration of power constraint at the relay, where 
= ^Pf^/ [(1 - P^){Ps\hs^? + + <j^) + 

The comparisons of secrecy rate performance achieved by 
PA and SAF are evaluated in Fig.2 with fixed Pd = 40dBm. 
For the special case W = 1, it is observed that SAF achieves 
a better performance than PA. The reason is that SAF can get 
the global optimal PS ratio by exhaustive search in this case, 
where PA just obtains a local optimal one, but the performance 
of PA approaches to be optimal in the whole region of SNRg, 
which indicates that PA is accurate and effective. When N = 
5, PA outperforms SAF in the whole region of SNRg. 


In Fig.3, we show the performances yielded by different 
number of relays and by different Pd resorting to PA and 
SAF with fixed Pg = 40dBm, respectively. When the power 
of AN is enhanced, both PA and SAF enjoy a better secrecy 
rate due to the weakened reception of eavesdropper. Moreover, 
it is noted that as the number of relays increases, more WPRs 
are participating in cooperative transmission yielding higher 
performances, meanwhile, PA achieves the maximum secrecy 
rate and the increasing performance gap with A^is notably seen. 

Fig.4 is presented to demonstrate the convergence result of 
PA, where we give the convergence result of PS ratio and the 
secrecy rate versus iteration times. Notice that the PS ratio p 
of each relay will finally converge to its own local optimal 
solution as the iteration times increases. 

V. Conclusions 

In this paper, we investigate the secrecy rate maximization 
problem in multiple non-regenerative wireless-powered relays 
networks based on power splitting scheme, which has not 
been studied yet in the literature. We propose an algorithm 
resorting to block-wise penalty function method to jointly 
optimize power splitting and beamforming for all active relays, 
which outperforms the benchmark. And we also prove that the 
proposed algorithm can converge to a local optimal solution. 
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